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ABSTRACT 

We present a continuing study of a sample 44 molecular outflows, observed in 
lines, closely associated with 6.7GHz methanol masers, hence called Methanol 
Maser Associated Outflows (MMAOs). We compare MMAO properties with those of 
outflows from other surveys in the literature. In general, MMAOs follow similar trends, 
but show a deficit in number at low masses and momenta, with a corresponding higher 
fraction at the high end of the distributions. A similar trend is seen for the dynamical 
timescales of MMAOs. We argue that the lack of relatively low mass and young flows in 
MMAOs is due to the inherent selection-bias in the sample, i.e. its direct association 
with 6.7GHz methanol masers. This implies that methanol masers must switch on 
after the onset of outflows (hence accretion), and not before a sufficient abundance of 
methanol is liberated from icy dust mantles. Gonsequently the average dynamical age 
of MMAOs is older than for the general population of molecular outflows. We propose 
an adjusted evolutionary sequence of outflow and maser occurrence in the hot core 
phase, where methanol masers turn on after the onset of the outflow phase. 

Key words: line: profiles; masers; molecular data; stars: massive, formation, proto¬ 
stars, outflows; submillimetre: stars 


1 INTRODUCTION 


The study of massive stars (> SM©) is of importance, as 
they serve as principal sources of heavy elements, UV ra¬ 


diation, and feedback in the interstellar medium (e.g. Zin 


necker & Yorke 20071. The formation and evolution pro¬ 


cess of young massive stars is complicated and not yet fully 
established. In order to improve our understanding of the 
important events which take place during the formation of 
high-mass stars it is necessary to identify and study different 
evolutionary phases. 


* Email: lientjiedv@gmail.com 


Our current understanding of the formation process of a 
massive star can be summarised into a crude four-stage evo¬ 
lutionary sequence (Zinnecker & Yorke 20071. Prior to the 
optically visible main-sequence life of OB-type stars, they 


spend 15% of their lifetime in an embedded phase (Church- 


well 20021. Zinnecker & Yorke (20071 sub-divide the em¬ 


bedded phase into four stages each of which is characterised 
by one or more types of objects. 


The study presented here focuses on the Hot Molecu¬ 
lar Core (HMC) stage, which is the first indirect manifesta¬ 
tion of a young massive star prior to the well-studied UCHll 
phase. HMCs are compact 0.1 pc), dense 10^cm“®), 
warm 100 K), and short-lived 10® yr) with high molec- 
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ular brightness temperatures (e.g. Garay & Lizano 1999 
Kurtz et al.[ 2000| Viti[ 20051. HMCs are warmed by an 
embedded (proto)star, but since the infalling material ab¬ 


sorbs the UV flux from this central protostar (Churchwell 


2002), they show no or only weak radio free-free emission 


(Beuther & Shepherd 2005 Dunham et al. 20111. 


HMCs are also characterized by significant abundances 
of large, complex organic molecules, including methanol (e.g. 
van Dishoeck fc Blake| 1998). During the initial cold (~ 10 
K) collapse phase, many molecular species that are present 
in the gas phase freeze out onto dust grain surfaces in the 
high-density inner region. These grains provide a site for 
further chemical evolution, producing more complex species 


(e.g. van Dishoeck & Blake 1998 Lintott et al.| 2005). It is 


believed that the complex chemistry in HMCs comes from 
these icy mantles that build up on the dust grains during 


the early collapse phase (e.g. Brown et al. 1988 Nomura & 


Millar 2004 Garrod fc Herbstf 2006[). 

Following this process, the associated protostar evolves 
to quickly warm the gas and dust (to ~ 100 K), evaporating 
the icy mantles, and injecting newly formed (generally more 
complex; [Lintott et al.j [2005ll grain mantle material into 
the gas phase (e.g. van Dishoeck fc Blake||1998| Viti et al. 


2004b Nomura fc Millarj 2004 Garrod & Herbst 2006). 
According to jViti et al.j ( |2004aP , this liberation of complex 
molecules into the gas phase is also assisted by molecular 
outflows developing from the central object, an important 


signpost of the HMC phase (e.g. Shepherd & Churchwell 
1996a||Cesaroni et al.||1997||Hunter et al.||1997||Kurtz et al. 


2000). These outflows clear out cavities in the envelope, driv¬ 


ing astrochemical processes via strong radiative (and possi¬ 


bly shock) heating of the cavity walls (e.g. van Dishoeck & 
Blake|[l998| ). 


Hence, this is the stage when large amounts of methanol 


are released from the grain mantles (Garrod & Herbst 


2006). The period from the heating of low temperature 
grains to the desorption of ices is fairly brief, ~ 10"^ — 10® yr 


(e.g. Bernasconi fc Maederj 19961. Viti et al. 12004b i found 
that organic molecules can be used as “chemical clocks” and 
showed that large species, such as methanol, are good indi¬ 
cators of the more evolved hot core phase, as these strongly 
bound species are abundant in the gas phase only at late 
times in a hot core’s evolution. They showed an increase 
in the fractional abundance of methanol by six orders of 
magnitude at ~ 2 — 6 x 10"^ yr since warm-up started, for 
10 —25M0 stars (see also |Viti]|2005|[Garrod fc Herbst[|2006| 


Garrod et al. 2008). 


Together with molecular outflows, interstellar 6.7GHz 
methanol masers are an important signpost of the HMC 
phase, being relatively common, the second brightest of the 
Galactic masers (only surpassed by 22GHz water masers), 
and radiating at cm wavelengths, hence not affected by 


high opacity gas and dust regions (e.g 

. Menten| 

1991 

Walsh 

et al. 1119981 Minier et al. 

20011 2003 

Beuther et al.l 

2002b 

Walsh et al.l 120031 Coc 

ella et al.l 

20041 Ellingsen 

2006 

Purcell et al. 

2006 Caswell 

2013). According to models 

by Sobolev et al. (1997); Cragg et al. (2001); Sutton et al. 

(2001|| and jCragg et al.j (|2005P, the radiative pumping of 


(class H) methanol masers by an infrared radiation source, 
happens in warm (100 K < Tdust), dusty environments, 
with gas densities between 10® — 10®cm“® and a high frac¬ 
tional methanol abundance of > 10“^ (column densities of 


~ 10^® — 10^^ cm“^) is required for 6.7GHz maser emission. 
HMCs have a brief lifetime, las ting betwe en 2 x 10® and 
6 X 10^ years (Kurtz et al. 2000). Work by Hartquist et al. 
11995) and van der Tak et al. (2000) showed that it requires 


10“^ — 10® yr to obtain a high methanol abundance via evap¬ 
oration of the icy dust grains in HMCs. Given the required 
abundances for methanol masers and methanol abundance 
models, it is thus clear that methanol masers are likely to 
occur towards the end of the HMC phase. Studying molec¬ 
ular outflows toward 6.7GHz methanol masers can reveal 
information regarding different sub-stages within the HMC 
phase, contributing to our developing understanding of the 
evolutionary timeline of massive protostars. 


In de Villiers et al. (2014 hereafter Paper I), high veloc¬ 


ity wings were identified in 58 three dimensional position- 
position-velocity cubes, each simultaneously observing the 
^®CO and C^®0 (J = 3 — 2) transitions with the HARP 
instrument on the JCMT, all targeted toward 6.7GHz 
methanol masers. Usable maps existed for 55 of these spec¬ 
tra, and were utilised to calculate the physical properties 
of the outflows. Relations between outflows and their as¬ 
sociated clump and maser properties were investigated and 
found to support the theory that massive star formation fol¬ 
lows a scaled-up version of low-mass star formation. Here, 
we compare the distributions of the physical properties of 
44 Methanol Maser Associated Outflows (MMAOs) (see 
with other CO outflow surveys. Without redefining the 
whole evolutionary cycle for high mass star formation, the 
work presented here addresses the relationship between the 
6.7GHz methanol masers, and the outflows they are associ¬ 
ated with. 


2 SAMPLE DESCRIPTION AND CAVEATS 

An initial sample of 70 6.7GHz methanol maser positions, 
drawn from the Methanol Multibeam (MMB) catalogue 


(Green et al. 2009 Breen et al., 2014 in prep.), were ob¬ 


served with the heterodyne array receiver, HARP, on the 
James Clark Maxwell Telescope (JCMT) on the summit 
of Mauna Kea, Hawaii. Targets were mapped in ^®CO and 
C^®0 (J = 3 — 2) transitions, covering a range of Galactic 
longitudes between 20.0° and 34°. A total of 58 spectra and 
their maps, extracted at peak ^®CO clumps, were of suffl- 
cient quality to be analysed. The complete list of these tar¬ 
get coordinates, and their 58 associated ^®CO peak clump 
emission coordinates, is given in Paper I, together with a 
description of the data reduction and analysis. 

All data cubes exhibited high velocity features, and 
^®CO integrated intensity maps were created using the blue- 
and red velocity wings. Outflow maps were created for 55 
targets (published in Appendix B of Paper I), and their 
physical properties tabulated (Table 5 of Paper I). A subset 
of outflows were selected based on their having resolved kine¬ 
matic distances and peak ^®CO clump positions within 18” 
(c.f. the 14” JCMT beam at these frequencies) from their 
associated maser positions. A total of 44 outflows satisfied 
these criteria, and were thenceforth referred to as Methanol 
Maser Associated Outflows (MMAOs). 

While it is not possible to resolve the contribution from 
individual stars or protostellar cores in our data, a relation¬ 
ship between the clump mass, the most massive star present 
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in the clump, and the associated outflow’s energetics can 
be assumed. As discussed Paper I, the higher mass accre¬ 


tion rates of the most massive and dense clumps | McKee & 


Tan 


20031 imply that they will arrive at the main-sequence 


ahead of their lower-mass neighbours (Mottram et al. 2011 


Urquhart et al. 20141. Thus, while the most massive star, 


traced by the 6.7GHz methanol maser in this study, is close 
to joining the main sequence, it is likely that the lower- 
mass members of the proto-cluster still need to evolve to a 
stage where they can significantly contribute to the observed 
clump luminosity. Results by Urquhart et al. (2013b I sug¬ 
gest that it is likely that the clump’s bolometric luminosity 
is actually dominated by the most massive stars (see Paper 
I for more detail). [Urquhart et al. (20141 also found that 
there is a strong correlation between the clump masses and 
the bolometric luminosities of the most massive stars, but 
that the total clump luminosities are much lower than would 
be expected from the fully formed cluster. This supports the 
hypothesis that the most massive stars have very rapid evo¬ 
lution times and are likely to dominate the observed clump 
luminosities, making it a reasonable assumption that so too, 
are the outflow luminosity and energetics. 


2.1 Outflow properties 

We compare the distributions of the outflow masses Mout, 
momenta p, and dynamical timescales ta with previous CO 
studies of massive molecular outflows, summarized in Ta¬ 
ble Although these samples are all biased in different 
manners, as described by their selection criteria, none of 
them are biased in terms of a 100% association with 6.7GHz 
methanol masers as is the case for MMAOs. One could ar¬ 
gue that the MMAO sample is the best defined, and most 
homogeneous of the four samples compared. However, we 
take note that the methanol maser selection was carried out 
towards maser sources with accurate positions at the time, 
biasing them towards the more luminous sources, detected 
prior to the MMB surveys in targeted searches towards e.g. 
IRAS selected targets ( jPestalozzi et al.||2005| and references 
therein). This brightness bias will be discussed in more de¬ 
tail in )]3.1| 

Figures to show the outflow mass, momentum and 
dynamical timescale distributions of MMAOs and the com¬ 
parison data. The calculation of these parameters are dis¬ 
cussed in detail in Paper I, and are briefly reviewed below. 
For fiv^ of the 6.7GHz methanol masers serving as trac¬ 
ers for MMAOs, trigonometric parallaxes have recently been 
published by jReid et al.| ( [^14| ). These distances are updated 
from the values given in Paper I (calculated using the Galac¬ 
tic Rotation Curve), together with their corresponding out¬ 
flow properties. These changes had an insignificant effect on 
the results presented in Paper I. 

The outflow mass is calculated as 

Mout = ^ (W X Ai) rriHa , (1) 

b,r 

where b and r represent the blue and red lobes respectively 
Ni the H 2 column density of each lobe, calculated from the 


1 G 23.437-0.184; G 23.706-0.198; G 29.865-0.043; G 29.956-0.016 
and G 31.282-1-0.062. 



log(Outflow mass / M©) 


Figure 1. Comparison between the distributions of the loga¬ 
rithm of outflow masses for MMAOs (grey shading) and other 
surveys (horizontal dash ed lines) (a) [Beuther et ah] l |2002a| | (b) 
|Wu et ar] ( |2004l l and (c) [Zhang et al.f < |2005| |. 


emission, with Ttrans = 31.8 K, the upper level energy 
of the J = 3 — 2 transition of ^^CO, and assuming an exci¬ 
tation temperature of Tex = 35 K (see Paper I). This 
column density is then converted to an H 2 column density, 
assuming a Galactocentric dependant isotopic ratio (jWilsonj 


& Rood 


19941, and the [CO]/[H 2 ] = 10 abundance ratio 
by Frerking et al. ( 1982| ). Ai is the surface area of each lobe, 
and mH 2 is the mass of a hydrogen molecule. 

The momentum was calculated per velocity channel 
(width Au) for each spatial pixel in the defined outflow 
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lobe area, using the channel velocity relative to the systemic 
velocity {vi), and the gas mass {Mi) corresponding to the 
emission in that channel. This is summed over all velocity 
channels, and all pixels in each lobe area, giving: 


p = An ^ 




b^Vi 


+E 


E 


The dynamical timescale td is calculated as 


td ■ 


In 


{Avb + AVr) /2’ 


( 2 ) 


(3) 


where Imax is the maximum lobe length between red (Ir) 
and blue (h) as measured from the clump coordinate to the 
furthest radial distance of the lobe. Avf, and Avr are the 
velocity extents, measured from the peak velocity as defined 
by the C^®0 spectrum, to the maximum velocity along the 
blue or red line wing (on a 3 cr level above the noise). 

In the case of monopolar outflow detections the length and 
velocity of the available lobe is used. As described in Paper 
I, the values for Avb and Avr are scaled by a factor two, to 
account for the difference in velocity extent between 
(this work) and the wider profile used in the compar¬ 

ison studies. For more detail regarding these measurements, 
we refer the reader to Paper I. 

It is interesting to note that the mass distribution for 
MMAOs shows a relatively smaller fraction of low-mass out¬ 
flows and a larger fraction of high-mass outflows compared 


to the other studies (except for the few outliers from Wu 
et al. |2004 1. The momentum distributions are shown ir 
Figure where the distribution for MMAOs have a simi¬ 
lar fraction to the general outflow population for momenta 
between IO^Mq km s“^ and 10^'^® Mq km s“^, but, as for 
outflow masses, low outflow momenta are under-represented. 

Figure shows the distribution for the dynamic 
timescale. Both [Beuther et al. (2002a I and Zhang et al. 
(20051 argue that for high-mass outflows, the dynamical 


timescale may often be a good estimate of the source age 
because the dynamical timescale of the outflow corresponds 
well to the free-fall timescale of the associated core. The age 
range of MMAO outflows is consistent with the other work 
considered here, with the interesting difference being a lack 
of the youngest outflows at < 5 x 10'* yr. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 
log(Outflow momentum / M© km s’') 


Figure 2. Comparison between the distributions of the loga¬ 
rithm of the outflow momenta for MMAOs (grey shading) and 
other surveys (horizontal dashed lines) (a) Beuther et al.| | ['2002a[ | 
(b) |Wu et al.| | [2004| | and (c) [Zhang et al.| < 2005| l. 


































































































































Table 1 . Other CO outflow studies used for comparison with the MMAO sample’s masses, momenta and dynamical timescales. 


BEUTHER ET AL. 


WU ET AL. 


ZHANG ET AL. 


KIM & KURTZ 


MMAOs 


Outflow sample size 


Used 18 out of 24, which had 
resolved distances. 


168 for Mout; 
139 for p; 135 for 


35 


11 


44 


i^CO(J = 1 - 0 ) and 

i2co(J = 2-l) 


Molecular tracer 


2CO(J = 2-1) 


2CO(J = 2-1) 


2CO(J = 2-1) 


3C0(J = 3-2) 


Facility 


30 m IRAM 
at Pico Veleta 


Literature study 


12 m NRAO at Kitt Peak & 
10.4 m CSO 


12 m NRAO at Kitt 
Peak 


15m JCMT at Mauna 
Kea 


Resolution 


Vel: 0.1 km smoothed to 1 
km s~^ 

Spatial: 11^^ 

(0.07 - 0.56 pc) 


Mixed 


Vel(NRAO): 

1.02 km (run 1 ) &; 
0.5 km (run 2) 
Spatial(NRAO): 29'' 
Spatial(CSO): 30" 
(0.04 - 1.57 pc) 


Vel: 0.68 km s ^ 
Spatial: 27" 
(0.06 - 0.60 pc) 


Vel: 0.5 km s ^ 
Spatial: 6 " 


Sensitivity 


Int.time 5 min.; rms: Not 
given. 


Mixed 


NRAO: Int. time = 4 min 
(run 1 ) 

Sz 1 min (run 2 ); 
rms = 0.2 K/2 km 
CSO: rms = 0.03 K/1.5 
km s~^ 


rms = 0.2 K 


rms = 0.24 K 


Selection criteria 


A sample of high mass 
protostellar objects from the 
IRAS catalog, following | Wood] 
I & Churchwell|l|1989] criteria, 
refined b y|Ramesh & | 

I Sridharan| l [l997] for the 

identification of UCHii 
regions. 


Combined study of 
391 high velocity 
molecular outflows 
from literature 
published prior to 
February 2003. 
High-mass outflows 
satisfies -Lout > IO^Lq 
and Mout ^ 3 M 0 . 


A selection from the 
flux-limited sample of bright 
IRAS sources, originally 260, 
following requirements from 
|Palla et al.|||1991], 101 


Molinari et al. (1998 


towards which IMolinari et al.1 

| |1996] found NH 3 emission, 
and |Molinari et al. | | |1998] 
observed 67 of them for radio 
continuum with VLA. |^hang | 
I et al. | ( |2005] used these 67, 
with additional two in their 
search. 


did not study 34 of the 
sources from [Molinari I 
I et al. | |1996] , of which 
Kim & Kurtz selected 
11 and added one 
additional source. 


Outflows associated 
with 6.7GHz methanol 
masers, where the 
latter have resolved 
kinematic distances 
and are located within 
18" of the peak ^^CO 
emission. 


References 


.Beuther et al. (2002a ; 


Sridharan et al. i 2002] 


Wu et al. (2004 


Zhang et al. 


(2005 


Kim & Kurtz 


(2006 


de Villiers et al. 


(2014; 


Notes 


The ambiguous distance to 
three sources was resolved 
using [Roman-Duval et al. j 
(120091. 


01 
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0 5 10 15 20 25 30 35 40 45 50 


Dynamical time scale t (x 10"^ yr) 


Figure 3. Comparison between the fractional abundance distri¬ 
butions of outflow dynamical timescales for MMAOs (solid lines) 
and other surveys (horizontal dash ed lines) (a) [Beuther et ah] 
) |2002a^ (b) |Wu et ah] | |2004[ | and (c) [Zhang et al.f l |2005[ |. 


3 DISCUSSION 

Considering Figures to there seems to be an under¬ 
abundance of MMAOs with low mass, momentum and dy¬ 
namical timscales. The latter is of particular interest. 

The two main independent variables in the calculations 
of outflow mass and dynamical timescale, are the outflow 
lobe’s projected size (included via the outflow’s area A^/b 
in the mass), or Zmax in the dynamical timescale, and the 
average velocity extent (Auavg = Au(, -I- Aur)/2. An impor¬ 



Max flow length (pc) 


Figure 4. Distributions of the maximum outflow length for each 
MMAO (shaded area) and the physical lengths for the resolution 
limit of a 14” JCMT beam given the distance to each MMAO 
target (red thick dashed line). 


tant question to answer is whether the comparatively high 
lower limit and high mean of MMAOs’ outflow lengths, is a 
resolution effect, or does it reflect an actual property of the 
sample? The same question can be asked about the lack of 
low mass and young outflows in our MMAO sample. Since 
the study has a 100% outflow detection rate (Paper I), these 
effects cannot be due to non-detections. However, since the 
full width half maximum (FWHM) of the JCMT beam is 
~ 14”, outflow lobes with an angular size smaller than this 
will not be resolved, imposing a limit on the smallest outflow 
sizes measured. In order to determine whether the observed 
distribution is affected by the telescope’s resolution, we (a) 
deconvolved the outflow lobe diameters with a 14” beam 
and re-measured the lobe lengths (Paper I), and (b) used 
the angular resolution limit of the JCMT as a hypotheti¬ 
cal physical lobe size at each MMAO’s distance. The distri¬ 
butions of the maximum deconvolved lobe length for each 
MMAO outflow (shaded area) and the linear size of a JCMT 
beam at each MMAO distance (red dashed line), are shown 
in Figure]^ The fact that the MMAO distribution shows a 
deficit in short outflow lengths, beyond the distance limit 
defined by the telescope’s resolution, indicates that such a 
deficit is not instrumental in nature, but is real. That is, if 
smaller outflows did exist in this sample of MMAOs, they 
would have been detected. 

The maximum lobe length ranges from 0.3 pc to 3.7 
pc for MMAOs, with a mean of 1.2 pc (see Table 3 in Pa¬ 
per I). These lengths are of a similar order of magnitude 
to both the samples of Beuther et al. 1 2002a|) and Wu et al. 


( |2004| ) , although they cover a larger range and have a higher 
mean than the outflow lengths from Beuther et al. (2002aI, 
which cover a range of 0.1 — 1.8 pc, with a mean of 0.8 pc. 
The MMAO values compare well with those from jWu et al.j 
(20041, ranging from 0.04 — 3.98 pc, with a mean of 1.1 pc. 
The distributions of these lengths are shown in Figure 1^ 
which does not show the outflow lobe lengths for the [Zhang 


et al. (20051 and Kim & Kurtz (20061 samples, as no lengths 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Max flow length (pc) 


25 


10 


JCMT beam limit 
MMAO deconv. (Tot. 44) 


L.l" 




10 20 30 

Dynamical time scale t (x 10^ yr) 


Figure 5. Distributions of the maximum outflow length for 
MMAOs (shaded are a), [Beuther et al.| ( |2002a| blue horizontal 
lines), and|Wu et al.|||2004| red diagonal lines). 


Figure 6. Distributions of the dynamical timescales for MMAOs 
after beam deconvolution (shaded area), and if every target only 
had a minimum lobe size of the 14" JCMT beam (red dashed 
line). 


are published by Zhang et al. (2005 I, leaving only the eleven 
outflow lengths published by Kim & Kurtz (20061 for that 
combined sample, which would not be statistically signifi¬ 
cant. Beuther et al. (2002a I detected a higher fraction of 


small outflow sizes compared to the MMAOs, which poten¬ 
tially contributes to the bottom-heavy distributions of their 
outflow masses and dynamical timescales. 

The average values and distribution of the outflow ve¬ 
locities, Auavg are found to range between 5 — 31 km s~^, 
with a mean of 13 km s“^ for MMAOs. iBeuther et al.l 
12002a I has a similar range and mean, from 4 to 2 7 km s“ 
with a mean of 13 km s“^. 


In 


Wu et al. 


(20041, a wider 


range of 3 — 133 km s“ is seen, but again a similar mean of 
18 km s“^ and median of 13 km s“^ (we calculated the lat¬ 
ter as the Wu et al. sample is skewed by a few high velocity 
outliers). Note that these values depend on the sensitivity 
of th e observations to outflow wings, but as mentioned in 
(2.1 our ^^CO measurements are scaled to equivalent val¬ 


ues expected in ^^CO (see Paper I). The velocity distribu¬ 
tions were found to be similar for all samples, and the fact 
that the outflow mass and momentum distributions are sim¬ 
ilarly top-heavy for MMAOs, implies that the velocity does 
not have much of an effect towards changing the shape of 
the MMAOs’ dynamical timescale distribution. No data is 
available for Zhang et al. (20051 and Kim & Kurtz (20061. 

Figure shows the distribution of dynamical timescales 
for our MMAO sample after the data have been deconvolved 
(the latter had minimal impact on the distribution), as well 
as the dynamical timescales calculated for minimum lobe 
length the size of the JCMT beam, corresponding to Figure 
It shows that our MMAO sample, as observed from the 
JCMT, is not sensitive (i.e. is incomplete) for ages below 
~ 2.5 X 10"^ yr. MMAOs show a turn-over at ~ 5 x 10'* yr, 
implying that the lack of the youngest outflows is not due 
to a resolution limitation. 

We also investigated whether a distance bias could 
cause the observed turn-over at low dynamical timescales. 
This is a risk because one tends to observe only the bright¬ 



Figure 7. Jackknife test for a Malmquist bias in the dynamical 
timescale distribution of MMAOs. 


est sources at the greatest distances (Malmquist bias). To 
investigate whether our data is affected by such a bias, we 
performed a jackknife test, ranking our dataset according to 
distance, and dividing it in half to form two distributions 
(Figure]^. Both the “near” and “far” samples show a turn¬ 
over in dynamical timescales at ^ 5 x 10* yr, indicating that 
it is unlikely to be a distance-dependent effect. 

Finally, we consider the inherent selection bias in our 
data, as all our targets are associated with 6.7GHz methanol 
masers. If the differences in the outflow mass and timescale 
distributions for MMAOs are due to this selection bias, one 
might argue that this makes our sample’s properties distinct 
from those of others. 

Although not a simple process, we expect that for any 
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individual outflow, mass will increase with time as long as 
accretion continues to power the expanding flow. The frac¬ 
tional distribution of MMAO mass favours large values, com¬ 
pared to the other samples (Figure [^, supporting the sug¬ 
gestion that the average age of MMAOs is older than a sam¬ 
ple not biased toward 6.7GHz methanol maser emission. 

The methanol maser association of our sample implies 
the presence of complex organic molecules in the gas phase 
at high column density. Viti et al. (2004b I found that large 
organic molecules such as methanol are good indicators of 
the more evolved hot core phase. Both |Hartquist et al(] 
119951 and Viti et al. (2004b I state that methanol shows 
a rapid increase to a high enough fractional abundance for 
the maser to switch on around 10^ — 10® years after in¬ 
fall and consequent warm-up commenced. This implies that 
there is a time in the hot core evolution, before which masers 
do not exist due to an insufficient fractional abundance of 
methanol. The turn-over seen in the MMAO age distribution 
in Figure is consistent with these ages. 

The evolutionary sequence of masers as given by |Breen| 
|et al.| | [2010[ their Figure 6) suggests that OH masers usually 
trace a more evolved stage of star formation compared to 
other masers (e.g. Garay & Lizano 19991. We used the cat¬ 
alogs in Caswell (20131 and Qiao et al. (20141 to match our 
MMAOs with OH masers, and found an association within 
~ 3.5” for at least 50% of the MMAOs. This is consistent 
with the suggestion that they occur at a later time in the hot 
core phase, although, given Breen et al. {2010l’s diagram, 
some 6.7GHz masers are also expected to occur earlier than 
the OH maser stage, which could account for the other 50% 
of our sample. 

Hence, we conclude that this study’s association with 
methanol masers introduces an age-bias towards more 
evolved outflows, suggesting that masers switch on some¬ 
time after the outflow has started. This implies that outflows 
found to be associated with 6.7GHz masers, should generally 
be slightly more evolved than massive outflows not associ¬ 
ated with these masers, and that it is unlikely that any of 
the very youngest (and probably also the smallest) outflows 
would appear in our sample, given its bias. 


3.1 Maser brightness bias 


The area of the sky from which our sample of MMAOs hap¬ 
pen to be selected includes a region of very high star forma¬ 
tion activity in the Galaxy, (i.e. mini-starburst W43 and the 
intersection of the end of the bar and one of the spiral arms). 
It is important to understand whether the sample could have 
been biased towards higher maser luminosities by this, and if 
so, to what extent this could influence the results. We found 
that the maser peak luminosities for MMAOs (Breen et al. 
in prep.), cover a range of ~ 2.5 orders of magnitude, while 
the MMB 6.7GHz masers cover a luminosity range up to ~ 6 


orders of magnitude ( Breen et al.[ 2011[ in their Figure 7). 

Although we do not sample below a 6.7G Hz maser 
luminosity of 100 Jy kpc^ (calculated as in Urquhart 
et al.[ |2013a[ their Figure 14), and an investigation of the 
Tmaser/Afciump distribution shows a bias to higher 6.7GHz 
maser luminosities for the MMAOs, our sample spans the 
luminosity range in which the bulk of the ATLASGAL pop¬ 
ulation lies. 

We investigate whether a correlation exists between the 



Figure 8 . Outflow mass as a function of the associated 6.7GHz 
maser’s luminosity for all MMAOs with available 6.7GHz lumi¬ 
nosity data (Breen et al. in prep.). 





Figure 9. Outflow dynamical timescales as a function of the as¬ 
sociated 6.7GHz maser’s luminosity for all MMAOs with available 
6.7GHz luminosity data. 


maser luminosity and outflow mass and dynamical timescale 
in Figures and Only a weak power-law relation exists 
for the outflow mass, Mout = (1.3 x 10^)L5;7 ghzi with a 
Spearman Rank coefficient of = 0.44 {ts = 2.96, where 
the minimum value of ts for the correlation to be significant 
is = 2.02 - see Paper I for more detail). It is thus possi¬ 
ble that the MMAO sample may have a slight bias towards 
more massive outflows, albeit with a weak correlation and 
substantial scatter in the data. 

Within our sample, no significant correlation is found 
between tdyn and Ls.tghz, with a power-law ht of tdyn = 
(1.15 X 10®)Le;7GHz (^s = 0.08 and ts = 0.45), although 
we note that our sample only covers the brighter end of the 
MMB luminosity function. Breen et al. ( 2010| ) proposed that 
the 6.7GHz maser luminosity increases as the source evolves. 
On the other hand, [Urquhart et al.| ( |2013a[ ) suggested that 
the strength of the maser emission is dominated by the en¬ 
ergy output of the central source and not driven by source 
evolution. Either way, any maser luminosity bias should be 
normalized out in tdyn calculations, as is seen in Figure 
since for a given age, the more energetic sources will have 
larger sizes (Zmax), but their outflow velocities (umax) will 
also be faster. This implies that the lack of “younger” out- 
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flows in Figure [^should persist regardless of a maser lumi¬ 
nosity bias. 


3.2 A revision to the evolutionary scheme for the 
hot core phase 

undertook a survey, searching for 

molecular outflows toward UCHll regions, for which they al¬ 
ready had 6.7GHz methanol and 22.2 GHz water maser data 
available. Based on the results, they proposed a tentative 
evolutionary scheme for the formation process of massive 
stars that suggests maser emission occurs before the outflow 
is developed enough to be detected. Following this, both the 
maser and outflow is present and detectable, after which 
maser emission disappears while the outflow is still present. 
Finally, only the UCHll region without masers or outflows is 
present. This sequence can be divided in the following steps: 


Codella et al. 


(2004 


(1) the earliest phase: maser emission present, outflow not 
yet developed enough to be detected, 

(2) both the maser and outflow are present and de¬ 
tectable, 

(3) maser emission disappears while the outflow remains, 

(4) outflow switches off, only the UGHll region remains. 


They argued that the lack of outflows in phase (1) does 
not imply a real lack of outflows, but instead reflects the 
youthfulness of the source whose outflow is too small to be 
detected by single-dish observations. Contrary to this stage 
where [Codella et al.| ( |2004| ) do not detect an outflow associ¬ 
ated with every 6.7GHz methanol maser, we found a 100% 
outflow detection rate towards 6.7GHz methanol masers in 
Paper I, which could be due to the higher sensitivity of our 
survey (the average rms of their J = 2 — 1 spectra is 

0.5 K versus our 0.24 K - see Paper I). Codella et al. (20041 
also showed that, if the outflow ages obtained in various 
survey s (Shepherd & Churchwell 1996a|b Beuther et al. 


2002aI, lying between 10 — 10° yr, are assumed to be rep¬ 


resentative estimates of their observed sample, a nd if the 
expected lifetime of an UCHll region is ~ 10° yr (Wood & 


Churchwell 19891, their time-line implies the maser phase 


should end after a few 10 yr, as the masers are expected 
to turn off once the expanding UCHll region dissipated the 


methanol below critical masing levels, (e.g. Walsh et al. 


1998 20031. 


From our results presented in §2.1[ supported by |Viti| 


et al. (2004b I’s chemical evolution models, we propose an 


amendment to this scheme, in which the maser switches on 
after the outflow has started and off before it disappears: 


(1) the earliest phase: outflow develops and grows with 
no methanol maser present; the infall process heats the dust 
grains and releases methanol into the gas phase, 

(2) both methanol maser and outflow are present and de¬ 
tectable (sufficient column density for methanol maser to 
turn on), 

(3) methanol maser emission disappears while the outflow 
remains, 

(4) outflow switches off, only the UCHll region remains. 


Thus, we propose that 6.7GHz methanol masers seem to 
signpost a more advanced stage of the hot core phase, and 
not as early as sometimes argued (e.g. [Codella et ^|2004| . 
Note that we only propose a change to the beginning of the 



Figure 10. Distributions of the dynamical timescales of the 
outflows from the surveys used for comparison in this study for 
whom 6.7GHz methanol maser association data were available. 
The grey shaded area shows the outflows with methanol masers 
associated, and the red dashed lines show the outflows without 
methanol masers. 


evolutionary scheme, not the end, which we have assumed 
stays unchanged. 


3.3 Other considerations 


Although MMAOs represent a homogeneous sample of ob¬ 
jects (i.e. observed with the same telescope and instrument), 
as with almost any statistical sample, there will be outliers 
and exceptions. Furthermore, these stages in the evolution¬ 
ary sequence of a massive YSO are not always clear-cut, and 
it is very possible for different evolutionary stages to overlap 


(e.g. Beuther et al. 20071. One example is G 24.790-1-0.083 
which is in the MMAO sample, but shows associations with 


Class I and H methanol, water and OH masers (Moscadelli 


et al. 20071. OH masers are usually indicative of a more 


evolved stage in the evolution of a massive YSO (Caswell 


19971, and yet this source has been calculated to have a 


relatively “young” age of ~ 2 x 10'* yr (Beltran et al. 


20111, younger than those calculated in Paper I. Similarly, G 
23.010-0.411 has an age of 1.8 x 10* yr (Sanna et al. 20141, 


and is also associated with Class I and H methanol, water 


20091. 


and OH masers (Forster & Caswell 1999 Cyganowski et al. 


We investigated whether a similar trend to that of the 
MMAOs in Figure [^ exists for methanol maser associated 
outflows in the literature. Outflows from Wu et al. (20041 


were cross-matched with the Methanol Multibeam Survey 
(MMB, [Green et al.j [2009^ , and The Arecibo Methanol 
Maser Galactic Plane survey (AMGPS, Pandian et al. 20071 


catalogues, with matching radii the size of the JCMT beam. 
Unfortunately, these two catalogues contain maser positions 
that were observed in the Southern Hemisphere, while most 
outflow surveys were observed in the North, implying a 
small cross-over area. Where information regarding 6.7GHz 
methanol maser association were available for the IBeutherl 
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et al. (2002a I; Zhang et al. (20051 and Kim & Kurtz (20061 


samples, we added them to the above two datasets. We found 
27 outflows associated with masers, aud 33 without. 

The distribution of these outflows’ dynamical timescales 
is shown in Figure |10[ Compared to outflows without as¬ 
sociated methanol masers, no deficiency in outflows with 
methanol masers is seen at young ages 2.5 x 10"^ yr), 
as is found for the MMAOs. However, a KS-test has shown 
that the null hypothesis that these two datasets are drawn 
from the same sample, cannot be rejected, and the difference 
seen between them is consequently statistically insignificant. 
We recognise that Figure]^ shows the presence of outflows 
associated with masers at young ages, whereas they appear 
to be absent for MMAOs at similar ages. However, we have 
shown that should these young outflows be present in the 
MMAO sample, they would have been detected (Figure]^, 
yet they are not. Furthermore, the samples plotted in Fig¬ 
ure are very heterogeneous (different observation tech¬ 
niques and calculation methods), while the advantage of the 
MMAO study is that observations and analyses were carried 
out in a homogeneous manner. 

In the near future, we should be able to assess these 
issues with stronger statistical support by expanding the 
MMAO sample to a broader, more representative range of 
6.7GHz methanol maser luminosities and Galactic longi¬ 
tudes. 


4 SUMMARY AND CONCLUSIONS 


The distributions of outflow masses, momenta and dynami¬ 
cal time-scales of MMAOs have been investigated and com¬ 
pared to previous related studies. We found that for an out¬ 
flow mass range between ~ 25 — 100 M©, and outflow mo¬ 
menta between ~ 100 — 500 Mq km s“^, the MMAO distri¬ 
bution follows a trend similar to that found in other studies. 
However, outside of these ranges, we find a comparatively 
smaller fraction of low-mass and momenta outflows and a 
higher fraction of high-mass and momenta outflows. 

The distribution of dynamical timescales shows that 
young MMAOs are few in number with an excess of old 
outflows when compared to other outflow surveys. This in¬ 
dicates a dynamically older sample of outflows compared to 
samples not biased by their association to 6.7GHz methanol 
masers. Given this, it suggests that masers and outflows do 
not develop at the same epoch, but that masers switch on 
after the onset of the outflow. This is consistent with IVitil 
I’s chemical evolution models which found that 
compounds, like methanol, are good indicators 
of the more evolved hot core phase, and that it is unlikely for 
6.7GHz methanol masers to occur before an age of 4 x 10^ 


et al. 12004b 


large organic 


yr- 

Furthermore, a modification to the first part of the 
evolutionary sequence for the hot core phase, presented by 


Godella et al. (20041, may be required. Without redefin¬ 


ing the whole evolutionary cycle for high mass star forma¬ 
tion, this study better defines the relationship of 6.7GHz 
methanol masers with outflows, where the masers appear 
at a more advanced stage of the hot core evolution than 
previously thought. 
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